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ABSTRACT 


A technique is presented for computing the 
matched-conic solutions of interplanetary trajec- 
tories for use as reference trajectories in guidance 
and navigation simulations. Continuity in position, 
velocity, and time at the sphere of influence is 
insured. The chief advantage of the matching tech- 
nique is its capability of being adapted to several 
mission types. Examples are presented of the types 
of missions considered: the free-flyby mission, the 
powered-flyby mission, and the stopover mission 
with a parking orbit about the target planet. 
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MATCHED-CONIC SOLUTIONS TO ROUND-TRIP INTERPLANETARY 

TRAJECTORY PROBLEMS THAT INSURE STATE-VECTOR 

CONTINUITY AT ALL BOUNDARIES 

By Victor R. Bond 
Manned Spacecraft Center 


SUMMARY 


A technique is presented for computing the matched-conic solutions of interplan- 
etary trajectories for use as reference trajectories in guidance and navigation simula- 
tions. The matching is made, with continuity in position, velocity, and time at the 
sphere -of- influence boundaries insured. The technique, which satisfies inflight con- 
straints at the target planet, is extended to several types of round-trip planetary mis- 
sions. The types of missions considered are the free-flyby mission, the powered-flyby 
mission, and the stopover mission with a parking orbit about the target planet. An ex- 
ample is presented of each type of mission. 


INTRODUCTION 


Matched-conic, or analytic, solutions to interplanetary trajectory problems have 
been used successfully in mission studies and as first approximations to more precise 
solutions. In reference 1, the need was stated for continuous conic solutions to inter- 
planetary trajectory problems. The primary use of such continuous solutions is in the 
area of guidance and navigation error analyses. A technique was also presented in ref- 
erence 1 for generating interplanetary trajectories that are continuous at the sphere of 
influence; however, the development was confined to one-way or single-leg trajectories. 
In manned interplanetary applications, the return trajectory must also be computed; in 
general, the departure and return portions of the trajectory cannot be computed inde- 
pendently. The two trajectories are related by conditions at the target planet that can- 
not be specified arbitrarily. 

Other matched-conic techniques for solving the round-trip interplanetary trajec- 
tory problem have been used. The circumlunar matched-conic trajectory technique 
presented in reference 2 and repeated later in reference 3, is probably one of the ear- 
liest techniques used. Several types of round-trip trajectories are discussed in refer- 
ence 4. Reference 5 is a discussion of round-trip missions and a presentation of a 
method for determining if an additional velocity increment on either the departure or 
the return trajectory, or on both trajectories, will reduce the total velocity require- 
ments for the mission. 



In this report, some of the mathematical techniques that are required for 
matched-conic trajectories to meet specified boundary conditions at the departure, tar- 
get, and return planets are presented. Boundary conditions at the planets are discussed 
in references 1 and 3. A detailed discussion of a process for matching the heliocentric 
and planetocentric trajectories at the sphere-of -influence boundaries will also appear 
in this report. A different matching technique was presented in reference 1. None of 
references 1 to 5, with the exception of references 2 and 3 which are concerned with 
circumlunar trajectories, present a detailed description of the matching of the depar- 
ture and return legs within the target-planet sphere of influence, although the assump- 
tion might be made that the matching is implicitly considered, if not explained. The 
matching of the departure and return legs will generally involve variation of either or 
both flight times until certain constraints at the target planet have been satisfied and 
the matching has been completed . 

The precision with which the matching at and within the spheres of influence is 
made is determined by the usage that the results of the solution will have. For studies 
in which only minimum velocity requirements are desired, little or no matching is re- 
quired for the single -leg trajectory. To find the round-trip trajectory with minimum 
velocity requirements, the trip times that satisfy mission constraints at the target 
planet must be provided first. The constraints arise from the type of mission (e. g. , 
free flyby, powered flyby, stopover with parking orbit, etc.). (The term ’ 'free flyby” 
refers to a trajectory in which no velocity change is applied in the vicinity of a planet 
while flying by the planet; ’’powered flyby" refers to a trajectory in which a velocity 
change is applied in the vicinity of a planet while flying by the planet. ) 

There are applications, however, for which a more precise matching is required. 
For guidance and navigation error analyses, it is desirable to have a reference or nom- 
inal trajectory that is continuous in both position and velocity. The matching technique 
described in this report is presented for use in guidance and navigation error analyses 
that require continuity in position and velocity. The equations of motion in guidance and 
navigation analyses are linearized by series expansions about the reference trajectory. 
The first-order terms of the series expansions yield the sensitivity coefficients for 
propagating errors along the reference trajectory. The sensitivity coefficients are 
discussed and derived in reference 3 and will not be considered in this report. 

In the section of this paper entitled "Example Problems, ” an example problem 
will be presented for each of the three types of trajectory solutions. The first two ex- 
amples were chosen from recently published literature. The author wishes to acknowl- 
edge E. W. Henry of the Manned Spacecraft Center for suggesting the third example 
problem. 


SYMBOLS 


a semimajor axis 

a.. matrix elements 

i] 

e eccentricity 
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f true anomaly 

f true anomaly of hyperbolic asymptote 

H argument of hyperbolic sine or cosine (analogous to eccentric anomaly) 

h altitude above planet surface 

A 

h unit vector along angular momentum 

i inclination with respect to planet equator (which may be arbitrarily defined) 

A A A 

i, j, k orthogonal - set of unit vectors 

J g coefficient of second zonal harmonic of planetary potential function 

k index specifying maximum or minimum periapsis declination 

m integer 

A 

n unit vector along ascending node of orbit and planet equator 

P computed period 

p semilatus rectum 

R d equatorial radius of a planet 

-D 

R position vector relative to Sun 

r position relative to planet 

r position vector relative to planet 

n 

A 

S unit vector along hyperbolic asymptote 

T Julian or calendar date 

t time from zero 

t stay time at the target planet 

s 

V velocity vector relative to Sun 

v velocity relative to planet 
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hyperbolic excess velocity 


v velocity vector relative to planet 

v oo hyperbolic excess velocity vector 

X, Y, Z Cartesian coordinates of the position vector 
y constraint function 

A 

a right ascension of S 

At specified time increment 

A 

5 declination of S 

St computed time increment 

i cos_1 (|) 


-1 

K = COS 


~ V °°,A,T 
i V °°» A, T 


V °°,D, t \ 


/i gravitational constant 

v one-half the angle between the arrival and the departure hyperbolic excess 

velocities at the target planet 



t designates either t R or t T in equations (54) to (59) 


< p angle between the vectors n and S 

fi right ascension of the ascending node 

u> argument of periapsis 


Subscripts: 

A arrival planet 

D departure planet 
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k index to resolve nodal ambiguity 

P planet 

R return planet (usually the same as the departure planet, but not required to be 

the same) 

T target planet 

it periapsis (or minimum distance from planet) 


Superscripts: 

* sphere -of-influence boundary 

' evaluated from heliocentric conic 

+ evaluated from planetocentric conic 

c computed quantity 

(j) quantity evaluated during jth iteration, where j = 0, 1, 

Operator: 

(’ ) derivative with respect to time 


COMPUTATIONAL TECHNIQUE 


The techniques discussed in this report have been programed in FORTRAN V for 
the Univac 1108 digital computer. Several types of trajectory problems are considered. 
The first problem considered will be to obtain the single-leg matched-conic trajectory 
between any two arbitrary planets. A flow chart describing the matching process for 
this trajectory solution is shown in figure 1. The single-leg matched-conic solution is 
required in order to obtain the round-trip trajectory solutions that will be discussed in 
the following order: (1) the free-flyby trajectory; (2) the parking-orbit (stopover) tra- 

jectory about an oblate planet; and (3) the coincident periapsis trajectory, which will 
yield either the powered -flyby trajectory or the parking orbit about a spherical homo- 
geneous planet. 

The computer program that has been developed has a limited search capability. 
That is, data such as departure date and approximate flight time are normally input to 
the program, and the input flight times are then changed by iteration until the con- 
straints, boundary conditions, and matching criteria are satisfied. The iterations 
usually occur in two modes, as shown in figure 2. TJie first mode, the gross iteration 
mode, utilizes heliocentric conics to satisfy constraints at the target planet. During 
the gross iteration mode, no matching of position and velocity vectors at the spheres 
of influence is achieved. After the constraints have been satisfied in the gross itera- 
tion mode, the tolerances are reduced, and the fine iteration mode begins. In the fine 
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Figure 1. - Computation for a single-leg 
(one way) interplanetary matched-conic 
trajectory. 



Figure 2. - Computation for round-trip 
interplanetary matched- conic trajec- 
tories. 


interation mode, single- leg matched conics are used to satisfy the same set of con- 
straints as were satisfied in the gross iteration mode. The end result of the fine iter 
ation is that the constraints and all the other boundary conditions are satisfied. 


DETERMINATION OF THE MATCHED-CONIC TRAJECTORY 
BETWEEN TWO ARBITRARILY CHOSEN PLANETS 


To obtain the matched-conic trajectory between two arbitrarily chosen planets 
(referred to as the single-leg matched-conic solution), T D is designated as the peri- 

apsis departure date, T A is designated as the periapsis arrival date, i D is desig- 
nated as the inclination of the trajectory with respect to a planeto centric coordinate 
system at the departure planet, i A is designated as the inclination of the trajectory 

with respect to a planetocentric coordinate system at the arrival planet, h^ D is des- 
ignated as the periapsis altitude at the departure planet, and h^ A is designated as the 

periapsis altitude at the arrival planet. To resolve the ambiguity of the trajectory in 
each planetocentric sphere of influence, kp is designated as an index that specifies the 

maximum (k = 1) or minimum (k = 0) periapsis declination at the departure planet, and 
k. is designated as an index that specifies the maximum (k = 1) or minimum (k = 0) 

dkV 

periapsis declination at the arrival planet. For a given set of the specified quantities 
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Tyv T . > i^, i*j h _ and h A , there are four possible solutions corresponding 
to the four possible combinations of and k^. The bounds on the specified inclina- 
tions i_^ and i. and the indices k^. and k. will be discussed in the section of this 
DA DA 

report entitled "The Planetocentric Phase." 


The Heliocentric Phase 

If the date of arrival or departure from a planet is given, then the heliocentric 
position and velocity of the planet may be computed from the planet ephemeris. At the 

departure planet, the position and velocity vectors are Rp and Vp ^T^j; 

and at the arrival planet, the position and velocity vectors are Rp ^T^ and 

Vp a(T^J. If the spacecraft is at a known position r^T^ or r^T^J from the 

departure or arrival planets at the same date, 
spacecraft are 


R d( T d) r D ( T d) 


and 


r a( t a) = ''a( t a) 


The position vectors R^ and R^ and the time T^ - T^ are used to determine 

the heliocentric velocities and in the vicinity of the departure and arrival 

planets. This calculation is known as Lambert’s problem and is documented in refer- 
ence 3. The velocities of the spacecraft with respect to the planet are then 


then the heliocentric positions of the 
+ R P, d( t d) ^ 

+ r p,a( t a) < 2 > 


V D* ( T d) ^d( T d) " ^P, D ( T d) 

and 

V( t a) = ? a( t a)- ? p,a( t a) 
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When the procedure is used for the first time in the solution of a single-leg tra- 
jectory, r^ and r^ are zero, and the spacecraft is assumed to be at the center of 
the planets at the dates and T^. After the planetocentric computations, good 

estimates for r D and r^ are obtained. The superscript * indicates that r^ 

and r^ are taken at the planet sphere of influence. 

The dates at the spheres of influence must also be corrected at arrival and de- 
parture. Therefore, the vectors 

«D ( T Dj= r D ( T D) +R P,D( T D) < 5 > 

and 

R A ( T Aj= r A ( T a) +R P,a( T a) < 6 > 

jjf jjc 

where T^ = T^ + t ^ and = T^ + t ^ are now used in the solution of 

— . + ak 

Lambert's problem to obtain V D and . The velocities at the spheres of influ- 
ence are 


and 







( 7 ) 


( 8 ) 


The superscript ' indicates velocities computed from heliocentric orbits, with the 
planets assumed to be massless. 
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The Planetocentric Phase 


To compute the trajectory within a planetocentric sphere of influence, the follow- 
ing vectors are specified: 

— * 

1. The position vector r of the spacecraft with respect to the planet at the 
sphere of influence was computed during the last sphere-of-influence computation and 
remained unchanged during the heliocentric phase. 

— * t 

2. The velocity vector v of the spacecraft with respect to the planet at the 
sphere of influence was computed from equations (7) and (8) during the heliocentric 
phase. 

3. The inclination of the hyperbola with respect to the planet equator (or to some 
other arbitrary plane) is designated as i. 

4. The periapsis altitude of the hyperbola is designated as h^. 

5. An index that specifies which node is to be chosen is designated by k. 

The computations for departure and arrival are similar, and the subscripts D 
and A are omitted, except where they are necessary for clarity. During the first 
computation of the trajectory in a planetocentric phase, only the magnitude of the posi- 
_** — * * 

tion r is known, and the velocity v must be approximated by the velocity of the 
spacecraft relative to the planet and must be evaluated at the center of the planet as- 
sumed massless, using equations (3) and (4). This evaluation causes no special prob- 
lem, since v^’ and v^* are good approximations to the sphere-of-influence 

5 jc f , j|c t sic sjc 1 

velocities v^ and v^ . The purpose for requiring both r and v is to 

allow the v^. vector to be computed. During the first computation in a sphere of 

influence, v^ must be approximated, although during subsequent passes, v rjo is 
computed exactly. 


When the spacecraft is assumed to be at the center of the planet ^ r D = 0 or 
r^ = o), an approximate computation of the hyperbolic asymptote is made from 


A 


S = 



(9) 


where the plus sign is used for departure and the minus sign is used for arrival. 
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The magnitude of the hyperbolic excess velocity is 



When the spacecraft is at the sphere of influence, the hyperbolic excess velocity 
and asymptote are computed exactly (ref. 3) from 



f^ - f >0 for departure, and f^ - f <0 for arrival. 
The hyperbolic excess velocity is 


v»=|'U (13) 

and the hyperbolic asymptote is 


A 

S = + 



( 14 ) 
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The trajectory within the sphere of influence is given initially by r and v and 
has orbital parameters that are not related to the specified orbital parameters. When 
the trajectory is propagated to infinity, using equation (11), the inclination and periap- 

m ^ afc f 

sis radius that would be found from r and v lose all meaning. Therefore, the 
specified periapsis conditions are now imposed. 

The only problem remaining in the sphere-of-influence phase is to determine the 
position and velocity vectors r^ and v^ at periapsis and the time t from periapsis 

A * 

to the sphere of influence, using the quantities v w , S, i, h^, and r . The vectors 

r^ and v^ and the time t^ are easily derived from figures 3 and 4, and the results 
for departure are stated as 


\ D - D^ OS - Sln ”( A X S d) 


(15) 


and 


V tt,D - v 7T, D 


sin T) S D + cos rj 


x s) 


(16) 




(a) Spherical geometry. 


(b) Planar geometry (h extending 
out of the page). 


Figure 3. - Relation between the periapsis vectors and the hyperbolic asymptotes 

at the departure planet. 
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For arrival, the results are stated as 


and 


where 


: »,A = r I r,A^ os, >®A +sin ^ X ®A) 


V A = V A 
7T ? A 7Tj A 


-sin rj + cos 7j 


( ! «a) 


r = h + planet radius 

7 T TT * 


% vH r ,r ' a . 




1 " 


COS 7 ] = 

e 


2 < V < 7T 


A A 


h = i sin O sin i - j cos sin i + k cos i 


(17) 


(18) 




> (19) 


J 


k 



(a) Spherical geometry. (b) Planar geometry (h extending 

into the page). 


Figure 4. - Relation between the periapsis vectors and the hyperbolic asymptotes 

at the arrival planet. 


12 


! 


The right ascension of the ascending node is ambiguous. From figures 3 

and 4 


! a - a 
a + v + v 


( 20 ) 


where a is the right ascension of the asymptote. For departure, the nodes that give 
the maximum (k = 1) and minimum (k = 0) periapsis declinations are 


— Ot +0 + 77 



( 21 ) 


and for arrival, the nodes that give the maximum and minimum periapsis declinations 
are 


= a - ct 


! 


— 01 + O’ + 71 


( 22 ) 


where 


. -1 
cr = sm 



(23) 


and 6 is the declination of the asymptote. Equation (23) indicates that the permissible 
inclinations are bounded by 7T-|6|>i>|6|. 

The time from periapsis to the sphere of influence is 


+ 


1 


3 

— (e sinh H - H) 


(24) 
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where the plus sign is used for departure, and the minus sign is used for arrival. Also 


H = cosh 


-1 



(25) 


The sphere-of-influence position and velocity vectors that satisfy the specified 
inclination and periapsis altitude are now computed (ref. 3) from 


and 


— * 
r 


+ */ 77 - 77 

= r — cos f + — 
r v 

7 7 7 7 


sin 


(26) 


where 


— * 
v 



-f- sin f + 

7 7 


(e + cos f) 

77 


(27) 


p = a(l - e 2 ) 

77 

2 < f < for departure 
77 < f < for arrival 


cos f = 


(^■ 1 ) 


»5 


(28) 


In equation (26), r has an additional superscript + to indicate that the position 

vector in equation (26) differs from the original r of equation (11). The vector r 
is used in equations (5) and (6) to determine the heliocentric position vectors of the 
spacecraft at the spheres of influence. 
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The Matching Process to Obtain a Single- Leg Trajectory 


The process of matching the position and velocity vectors at the spheres of influ- 
ence to obtain the single- leg matched conic is accomplished by using a successive ap- 
proximation technique. The logic of the computation technique is also presented in 
figure 1 and in the following: 

1. If the departure and arrival dates T D and T^ are given, the relative veloc- 
ity vectors v^’ and v^’ are provided by the heliocentric solution, with the space- 
craft assumed to be at the center of the departure and arrival planets at the dates T^ 
and T^, respectively. 

2. If the relative velocity vectors v^' and v^' are given or if the sphere-of- 

^ * f * t 

influence velocities v^ and v^ are computed in the heliocentric phase, along 
with the inclinations i^ and i^ and the altitudes h^ ^ and h^ the position and 

— * _L * 

velocity vectors r^ ’ v d ’ r A ’ anc * V A s P^ eres °f influence are com- 

puted in the planetocentric phases. The times t^ ^ and t^ ^ from each periapsis to 
each sphere of influence are also computed. 

3. The heliocentric phase is now repeated, adjusting the heliocentric positions 
of the spacecraft and the departure and arrival times according to equations (5) and (6). 

* T * t 

The velocities v^ and v^ at the arrival and departure spheres of influence, 
respectively, are obtained from the adjusted solution. 

I _ jjc f ^ | I _ ^ j|e f A 

Vjj “ v d V A “ V A are now C0m puted. 

If the errors are less than the tolerance, the solution is assumed to be converged. If 
the tolerance is not met, then steps 2 and 3 are repeated. 

ije 

The matching process may be visualized by considering that the date T and the 

position vector r at the sphere of influence are changed until v = v . This pro- 
cedure for matching is similar to the procedure used in reference 1. 


THE FREE- FLYBY MATCHED-CONIC MODE 


To compute the orbital parameters for the free-flyby mode, is designated 

as the periapsis departure date, i^ is designated as the inclination at the departure 
planet, i^ is designated as the inclination at the return planet, h^ ^ is designated 
as periapsis altitude at the departure planet, h „ is designated as periapsis altitude 

7 T. rt 
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at the return planet, h T is designated as the periapsis altitude at the target planet, 

kp is designated as the maximum (k = 1) or minimum (k = 0) periapsis declination at 

the departure planet, and k R is designated as the maximum or minimum periapsis 

declination at the return planet. The flight times t T (time from departure to target) 

and t_ (time from target to return) are dependent variables and require initial guesses. 

An additional constraint must be satisfied in order to attain the free-flyby trajectory: 
the velocity magnitudes of the arrival and departure hyperbolic excess velocities must 
be equal. The flight times must be computed by a numerical iteration procedure so 
that the flyby constraint is satisfied and the computed periapsis altitude is equal to the 
specified value. 

In the first iteration, only heliocentric conics are used, and all the boundary con- 
ditions except Tp and h^ ,j, are ignored. This mode is called the gross iteration 

mode and is done to improve the flight times, to determine the flyby inclination hp, 

and to resolve the nodal ambiguity. 

When the flyby constraints are approximately satisfied, the computation proceeds 
to the fine iteration mode. In this mode, all boundary conditions, flyby constraints, 
and continuity at the spheres of influence are satisfied. The trajectories used in this 
mode are single-leg matched-conic trajectories, the solutions of which are described 
in a previous section. The computation for both the gross and the fine iteration modes 
is illustrated schematically in figure 2. 


Computation of the Free-Flyby Trajectory Within the 
Target- Planet Sphere of Influence 


It is assumed that the hyperbolic excess velocities v^ ^ and v^ ^ rp have 

been computed and that the magnitudes of these velocities are approximately equal. The 
hyperbolic excess velocities may be computed by either of the methods described in the 
section of this report entitled M The Planeto centric Phase, *' depending upon whether the 
velocities are relative velocities and the planets are assumed to be massless or upon 
whether the position and velocity vectors at the sphere of influence are known. 

As shown in figure 5, the plane of the orbit is computed from 


£ _ V °°, A, T X V °°, D, T 

T " - 

V «>, A, T V °°, D, T 


(29) 
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Figure 5. - Geometry of the free-flyby 
hyperbolic trajectory. 


The inclination is found from 


cos i T = & • A t 


( 30 ) 


The angle 2v between ^ ,j, and 

v _ _ is found from 
D, T 


sin 2 v = 


V< *>, A, T X v °°, D, T 
( V °°, A, t)( V °°, D, t) 


(31) 


The radius of periapsis is computed from 


T = a(l - esc v) (32) 


where 


a = 


M - = _ Ul 
( V °°) D, t) ( V °° 5 A, t) 


(33) 


It should be emphasized that r _ c is a computed value because in the free-flyby 

TT, 1 

mode, r T is actually specified. In the description of the iteration technique given 

TT, 1 

in the following section of this report, details will be given concerning how the flight 
times t T and t R are adjusted so that the difference r 7r ^ ~ r 7T T c vanishes. 

The inclination i T is to be used in the single-leg computations in place of a 
specified i,p. An additional quantity k is required in order to determine which am- 
biguous node is to be chosen. The decision about which k to use is resolved by ex- 
amining figures 3 and 4. 
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The nodal vector is 


, ix6 T 

n. = — : 

k sin i T 


(34) 


For the departure leg 


c°s = n t 


D,T 


(35) 


and cos must always be negative for a maximum (k = 1) periapsis declination and 
positive for a minimum (k = 0) periapsis declination. If cos <0, = a + o + t, 

if cos > 0, = a - or. Similarly, for the arrival leg 


cos $ A 




(36) 


The cos <£> A must always be positive for a maximum (k = 1) periapsis declination and 
negative for a minimum (k = 0) periapsis declination. If cos 0 A > 0, = a - <j; if 

cos 0 A < 0, = a + a + 7T. 


Iteration Technique for Computing the Flight Times 
for the Free-Flyby Trajectory 

In the solution for the free-flyby trajectory the flight times t T and t R must 
satisfy the constraints 


y l " V «,A,T " V °°,D,T ° 

(37) 

y 2 = h 77, T " h 7T, T = 0 

(38) 
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The periapsis altitude h „ C is computed from equation (32), and h is specified. 

7T ? ± " ) 

During the gross iteration mode, when only the constraints in equations (37) and (38) 
are to be satisfied, velocities that reduce the final errors in equations (37) and (38) to 
0. 1 km/hr and 0. 1 kilometer, respectively, are provided from equations (3) and (4). 
When this accuracy has been obtained, the fine iteration begins, and the flyby con- 
straints and all other boundary values must now be satisfied. In the fine iteration 
mode, the more precise velocities computed from equation (10) are used. The fine 
iteration mode is continued until the errors in equations (37) and (38) are 0. 01 km/hr 
and 0. 01 kilometer, respectively. 

The variables yj and y ^ may be considered functions of the flight times t T 
and tj^. 


*1 



(39) 


*2 



(40) 


The solution to these equations may be done numerically by the Newton- Raphson 
method, as described in reference 6. Symbolically, the solution may be written as 


t 


T 



+ §t T 


(41) 


and 




(0) 


+ 6t 


R 


(42) 


where 


'6t„ 


6t. 


R/ 


a ll a 22 " a 12 a 21 


22 


-a 


21 



(43) 
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The matrix elements 


a.. 


i] 


are defined by the partial derivatives 


ii " at T 
ay 2 

a 2i " at T 


Hi 

*12 ” 

d j2 

a 22 “ at„ 

K 


\ ( 44 ) 


The derivatives are difficult to obtain in closed form, but they may be obtained numer- 
ically by the following procedure: 

1. A pair of trajectories (that is, one trajectory from the departure planet to 
the target planet and another trajectory from the target planet to the return planet) is 

computed, using the first-guess times t T ^ and t R ^°\ Then, the residuals 


*1 


( 0 ) 



(45) 


and 



(46) 


are formed. 

2. A second pair of trajectories is computed, using the times 
t ^ = t T ^ + At T and t R ^; that is, t T ^ is incremented by At^, and t R ^ 
is held constant. Then, the residuals 
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and 


y W-y 17 < X >t 

y 2 y 2 |_T ’ R _ 


( 48 ) 


are formed. 

3. A third pair of trajectories is computed, using the times t T ^ 

t_/^ = t_,^) + At,,, and the residuals 
rt rt rt 


y i 


( 2 ) 



and 


( 2 ) 


= y 2[lT 


( 0 ) . ( 2 ) 

R 


(49) 


(50) 


are formed. The partial derivatives may now be approximated. 


y (1) _ y (0) 

y l y i 


11 " At n 


(1) (0) 

y 2 " y 2 


21 “ At n 


y ( 2 } _ y (0) 

y l y l 


l 12 ~ At 


a 22 “ At 


R 

( 2 ) _ ( 0 ) 

y 2 


R 


•> 


> 


(51) 


4. A fourth pair of trajectories is computed, using the corrected times from 
equations (41) and (42). If equations (37) and (38) are still not satisfied to some toler- 
ance, then the entire procedure is repeated. 
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MATCHED- CONIC TRAJECTORY WITH A PARKING ORBIT 
ABOUT AN OBLATE PLANET 


To obtain a matched- conic trajectory with a parking orbit about an oblate planet, 
Tp is designated as the periapsis departure date, T R is designated as the periapsis 

return date, T T is designated as the periapsis target date, ip is designated as the 

inclination at the departure planet, i R is designated as the inclination at the return 

planet, h ff R is designated as the periapsis altitude at the departure planet, h^ R is 

designated as the periapsis altitude at the return planet, h is designated as the 

“9 I 

periapsis altitude at the target planet, t is designated as the stay time at the target 
planet, kp is designated as the maximum (k = 1) or minimum (k = 0) periapsis decli- 
nation at the departure planet, and k R is designated as the maximum or minimum 

periapsis declination at the return planet. In addition, the constraint that the stay time 
be equal to an integral number of orbital periods must be imposed. 

It is assumed that the hyperbolic excess velocities v co a T arU * d T ^ ave 
been computed. These velocity vectors were computed, with a specified stay time t 

s 

at the target planet assumed. If the radius of periapsis at the target r is also 

7T, I 

specified, then it is possible to find several parking-orbit solutions that contain 

v . „ initially and that contain v^ n at the end of the stay time (fig. 6). The 

rotation of the parking orbit is accomplished by the secular perturbations that arise 
from the oblateness of the target planet. From reference 7, the rates of change of 
the node and of the argument of periapsis of the parking orbit are 



Figure 6. - Geometry of the parking orbit 
at the time of arrival and departure. 
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where J ^ is the second zonal harmonic of the planetary potential function and Rg is 
the equatorial radius of the planet. 

Not only are the solutions nonunique; they must be found numerically. The tech- 
nique used to obtain the solutions is discussed at length in reference 8 and will not be 
repeated here. It is sufficient to say that the technique yields several pairs of inclina- 
tions and eccentricities for the parking orbit. 

The solutions obtained in reference 8 did not impose the constraint that the stay 
time be equal to an integral number of orbital periods. When this condition is imposed, 
a unique solution is obtained. The condition may be expressed as 


t - mP c = 0 (54) 

s 

/ c c 

P . The computed period P is the anomalistic 

period, that is, the time required to go from periapsis to periapsis. Since the mean 
motion also undergoes a secular perturbation, the anomalistic period is not the same 
as the period in the corresponding two -body problem. From reference 7, the perturbed 

mean motion 2v/ P c and thus the computed period are found from 



2a' 


3J 2 R B 

Kl - e 2 ) 


3/2 


/3 . 2. 
\2 SU1 1 


‘) 


(55) 


The constraint in equation (54) can be satisfied either by permitting t^ to be free and 
constraining the total flight time, or by permitting t^ to be free and constraining t^. 

If t is used to designate either t^ or t^,, the constraint given by equation (54) 
may then be written as 

y( T ) = t - mP C = 0 (56) 

s 
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The time r may be found iteratively by the Newton-Raphson technique 


T U)- T (0> y (0) (T) 

dy 
dT 


where 


dy ,W-,W 
dT At 


( 57 ) 


(58) 


The Newton-Raphson technique requires two pairs of trajectories in order to 
establish the derivative in equation (58). The first iteration is a gross iteration; that 
is, no attempt is made to solve all boundary conditions or to force continuity at the 
spheres of influence. Succeeding iterations are fine iterations requiring all boundary 
conditions to be matched and requiring continuity at the sphere of influence. The only 
boundary condition that is not specified is the inclination at the target planet, which is 
found from equations (52) and (53). 


MATCHED-CONIC TRAJECTORY WITH COINCIDENT PERIAPSIS 
POSITIONS AT THE TARGET PLANET 


To obtain the matched-conic trajectory with coincident periapsis positions at the 
target planet, the same specifications are required as were required in the previous 
section. If the stay time is not equal to zero, then the solution yields the parking orbit 
about a spherical homogeneous planet, with the further restriction that the stay time, 
if not zero, must not be less than the stay time of a circular orbit at the specified alti- 
tude h If the stay time is specified to be zero, then the solution yields a powered- 

7T, 1 

flyby trajectory with an impulse applied at the periapsis. Both modes may be handled 
conveniently by the same information, since both modes have the common constraint 
that the periapsis position vectors for arrival and departure at the target be coincident. 
As in the previous section, either t,^ or t^ may be left free so that the periapsis 

altitude h „ may be specified. 

IT, X 

It is assumed that the hyperbolic excess velocity vectors v TO . T and v w n T 

have been computed. The inclination and the choice of the node to be used are deter- 
mined exactly as in the computation of the free-flyby trajectory. The radius of peri- 
apsis must be computed differently, since the matched-conic trajectory with coincident 
periapsis positions at the target planet is discontinuous in velocity at periapsis. 
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The solution for the radius of periapsis starts with the requirement that the peri- 
apsides of the arrival and departure hyperbolas must be coincident. In figure 7, it is 
shown that this requirement may be expressed as 



+ K = 2,r ( 59 ) 


/ A 



Figure 7. - Coincident periapsides 
condition. 


where 


-v 


K = COS 


°°, A, T °°,D, T 
v °°, A, T V °°, D, Ty 


(60) 





-1 


(Vt )( v °°,a,t) 


2 




(61) 


and 



For the given flight times t T and t R and stay time t g , the computed radius of 

periapsis r^ T c is a function of the velocity vectors v^ A T and T only. 

Equations (59) to (62) may be solved iteratively for r T C by using a Newton- Raphson 

"9 1 

technique. 

A second iteration is required in order to satisfy the constraint 


h 


rr,T 



(63) 


This constraint may be satisfied either by permitting t T to be free and constraining 
the total time, or by permitting t R to be free and constraining t T> The computation 
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is similar to the computation for the parking orbit about an oblate planet, which is dis- 
cussed in the previous section of this report. 


EXAMPLE PROBLEMS 


Three examples were chosen to illustrate the use of the matching technique dis- 
cussed in this report. In each example, the departure date and launch times were 
chosen from recent literature or from unpublished studies. The flight times were ad- 
justed during the computation process in order to satisfy the mission constraints. The 
purpose of the examples is not to present a comparison between the flight times chosen 
from the literature and the flight times computed by the iteration techniques discussed 
previously, but to demonstrate that continuous matched-conic solutions to round-trip 
trajectory problems can be generated from flight times that are only approximately 
known. 


The first example, presented in table I, is an Earth-Mars-Earth free-flyby tra- 
jectory. The data for this trajectory were chosen from reference 9. The initial 
guesses for flight times were found to be 130 days from Earth to Mars and 540 days 
for the return to Earth for an Earth departure date of Julian date 2 442 670. In table I, 
it is indicated that these times were changed to 136. 61225 and 539. 64008 days, respec- 
tively, in order to satisfy the constraints. The flyby and altitude constraints presented 


-3 -4 

in equations (37) and (38) were satisfied to 10 fps and 10 nautical miles, respec- 
tively. The boundary conditions were = 35 


*R = 50 ’ 


h _ = 262 nautical miles, 

n, D 

h „ = 50 nautical miles, and h — = 200 nautical miles. The periapsis at departure 

7 Ty I\ 71, 1 


and at return was chosen for a minimum declination 


(“o ' 0 


= 0 and k T 


= °> 


The second example, presented in table n, was chosen from reference 10. This 
stopover trajectory is of the low-energy, conjunction-class of Mars missions. The 
flight times were chosen to be 302 days from Earth to Mars and 321 days for the return 
to Earth. The Earth departure date is Julian date 2 444 930, and the stay time at Mars 
was chosen to be 416 days. The oblateness of Mars was used in order to rearrange the 
orientation of the parking orbit at its periapsis so that no plane changes were required. 
The time to return was allowed to be free in order to satisfy the integral orbit con- 
straint presented in equation (56). The return time was changed to 321. 8555 days, and 
the chosen parking-orbit solution had an eccentricity of 0. 781518 and an inclination to 
the Mars equator of 146. 580°. The boundary conditions were = 70°, = 75°, 

h ~ = 262 nautical miles, h n = 20 nautical miles, and h _ = 200 nautical miles. 

TT, JJ 71, rt 7T, X 

The periapsis at the departure planet was chosen to have a minimum (k^ = 0) declina- 
tion, and the periapsis at the return planet was chosen to have a maximum ^k^ = 1) 
declination. 


The third example, presented in table III, is an example of a parking orbit with 
coincident periapsis of the arrival and departure legs at the target planet. The depar- 
ture planet is Mars, the target planet is Venus, and the return planet is Earth. The 
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initial guesses for the flight times were 260 days from Mars to Venus and 170 days 
from Venus to Earth. The departure date at Mars was chosen to be Julian 
date 2 443 493, and the stay time at Venus was chosen to be 10 days. The total flight 
time was held constant, and the flight times were changed to 270. 68412 and 
169. 31588 days in order to satisfy the constraint, presented in equation (63), for the 
specified periapsis altitude. In order to obtain a highly elliptical orbit, the period of 
the parking orbit was chosen to be one-third of the specified stay time. The boundary 
conditions were i_ = 40°, i_ = 40°,. h ^ = 150 nautical miles, h = 50 nautical 

miles, h _ = 200 nautical miles, 1^ = 0, and k = 0. 


CONCLUDING REMARKS 


A technique for matching conic trajectories at the gravitation sphere -of -influence 
boundaries is presented. The matching is made, with continuity in position, velocity, 
and time at the sphere-of-influence boundaries insured. The technique is extended to 
several types of round-trip planetary missions and has the capability of satisfying in- 
flight constraints at the target planet. The types of missions considered are the free- 
flyby mission, the powered-flyby mission, and the stopover mission with a parking 
orbit about the target planet. An example is presented of each mission type. 

The chief advantage of the technique is its capability of being adapted to several 
mission types. The single-leg matched-conic trajectory can be used to obtain other 
types of solutions, for example, the optimum powered-flyby trajectory. The require- 
ment for its use in obtaining other solutions is that the constraints be properly identi- 
fied and stated mathematically. The limitations of this technique are no more severe 
than those for any other matched conic; that is, the technique represents only an ap- 
proximation to the precision-integrated trajectory. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, August 19, 1968 
981-30-10-00-72 
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TABLE I. - EARTH- MARS-EARTH FREE- FLYBY TRAJECTORY FOR AN EARTH DEPARTURE DATE 
OF SEPTEMBER 14, 1975 (JULIAN DATE 2 442 670. 0) a 


Location 

(coordinate 

Time from 
departure 
periapsis, 
days 


Position vector 
(b, c) 



Velocity vector, fps 
(c) 


system) 

X 

Y 

Z 

Magnitude 

X 

Y 

Z 

Magnitude 

Earth periapsis 
(planetocentric) 

0. 0000000 

2885.0346 

-1531.3456 

-1750. 8615 

3705.9361 

23741.259 

29758.012 

13093. 283 

40256. 933 

Exit Earth sphere 
of influence 
(planetocentric) 

1.7417041 

-7409. 3675 

417898.25 

273257.02 

499362. 84 

-580.81599 

16264.951 

10712.090 

19484.220 

Exit Earth sphere 
of influence 
(heliocentric) 

1.7417041 

0. 99728074 

-0. 12075733 

0. 00104548 

1. 0045657 

10155.239 

115785.87 

3357.2698 

116278.84 

Enter Mars 
sphere of 
influence 
(heliocentric) 

135.82706 

-0. 39312709 

1. 5504021 

0. 042567679 

1. 6000336 

-78956.047 

14540. 547 

59. 655663 

80283. 798 

Enter Mars 
sphere of 
influence 
(planetocentric) 

135.82706 

-303135.60 

-51517.335 

52282.088 

311895.24 

27037. 989 

4771.2949 

-4550.4197 

27830. 278 

Mars periapsis 
(planetocentric) 

136.61225 

160. 34837 

1736, 1971 

1066.6065 

2043.9524 

30935. 784 

1786. 1843 

-7558.2427 

31895.771 

Exit Mars sphere 
of influence 
(planetocentric) 

137.39745 

296818.96 

-16876.856 

-94299.085 

311895. 24 

26441.853 

-1683.4489 

-8515.7992 

27830. 278 

Exit Mar s sphere 
of influence 
(heliocentric) 

137.39745 

- 0.41436130 

1.5542281 

0.41691174 

1. 6090551 

-74638. 109 

12882,585 

-6356. 1945 

76007. 956 

Enter Earth 
sphere of 
influence 
(heliocentric) 

675. 14776 

0.47414165 

-0. 90459878 

-0, 00233758 

1.0213299 

83433.957 

74223.266 

11706.371 

112282.49 

Enter Earth 
sphere of 
influence 
(planetocentric) 

675. 14776 

27632.274 

-348254. 19 

-356817.50 

499362.84 

-1432.8183 

21953.590 

22279.237 

31310.980 

Earth periapsis 
(planetocentric) 

676.25233 

2790. 5805 

2099. 9400 

i 

102.48199 

3493.9361 

-19392.789 

23980.287 

36689. 108 

I 

| 47929. 375 


CO 

CD 


a For this trajectory, = 35°, = 50°, D = 262 n. mi., R = 50 n. mi. , T = 200 n. mi. , k^ = 0, and k R = 0. 

^Planetocentric position vectors are measured in nautical miles. Heliocentric position vectors are measured in astronomical units. 

c Planetocentric vectors are in the -planetocentric equatorial (Earth equatorial at Earth, Mars equatorial at Mars) system. Heliocentric vectors are in the heliocentric ecliptic 
system. 



TABLE n. - MARS STOPOVER TRAJECTORY FOR AN EARTH DEPARTURE DATE OF NOVEMBER 21, 1981 (JULIAN DATE 2 444 930. 0) 5 


1 

Location 

(coordinate 

Time from 
departure 
periapsis, 
days 


Position vector 
(b, c) 

— 


Velocity vector, fps 
(c) 


system) 

X 

Y 

Z 

Magnitude 

X 

Y 

Z 

Magnitude 

Earth periapsis 
(planetocentric) 

0 . 0000000 

2203.8515 

1089. 0369 

-2773. 2652 

3705.9361 

-29481. 566 

6540. 4417 

-20859.963 

36702.563 

Exit Earth sphere 
of influence 
(planetocentric) 

3.1134812 

-462994. 31 

-78452.616 

169837. 30 

499362.84 

-9444. 9839 

-1700. 8943 

3740. 5277 

10300. 112 

Exit Earth sphere 
of influence 
(heliocentric) 

3. 1134812 

.45 39 2 720 

. 87354679 

. 00231538 

.98444870 

-97533.951 

45080.633 

4108.2831 

107526.80 

Enter Mars sphere 
of influence 
(heliocentric) 

299.86212 

.18029548 

-1.4286007 

-.03421212 

1.4403391 

72089.638 

14845.003 

-2083.5210 

73631.725 

Enter Mars sphere 
of influence 
(planetocentric) 

299.86212 

64926. 100 

-270147.78 

141715.98 

311895.24 

-2190. 0579 

8626.0528 

-4586.8241 

10012. 197 

Mars periapsis 
(planetocentric) 

302. 00000 

-1729.9204 

320.21172 

-1040.4716 

2043.9524 

5549.2527 

17229.915 

-3923.7452 

18521.877 

In parking orbit 

302.00000 

-1729.9204 

320.21172 

-1040.4716 

2043. 9524 

4420.7970 

13726.165 

-3125.8409 

14755.403 

Out of parking orbit 

718.00000 

-626.75358 

1646.2809 

1036. 6680 

2043.9523 

14016.807 

3326.9259 

3191.0204 

14755.403 

Mars periapsis 
(planetocentric) 

718.00000 

-626.75358 

1646.2809 

1036.6680 

2043.9523 

16517.016 

3920. 3572 

3760.2100 

17387.357 

. Exit Mars sphere 
of influence 
(planetocentric) 

720.73968 

284451.20 

-115326.59 

-55370.858 

311895.24 

6988. 1476 

-2937.4740 

-1428.7832 

7713.9083 

. Exit Mars sphere 
of influence 
(heliocentric) 

720. 73968 

-1.4095229 

. 88520748 

.05195948 

1.6652467 

-36745.480 

-53721.482 

-2739.5946 

65143.943 

Enter Earth sphere 
of influence 
(heliocentric) 

1037.95365 

.99805264 

.02744371 

-.00191360 

.99843170 

7751.2668 

108673.38 

5767.1951 

109102.00 

Enter Earth sphere 
pf influence 
(planetocentric) 

1037.9537 

-350685.71 

-232226. 83 

-269171.84 

499362.84 

12405.649 

8111.4572 

9803.2147 

17770.731 

Earth periapsis 
(planetocentric) 

1039.8555 

979.42250 

-294.46639 

3309. 512f 

3463. 9362 

28086. 748 

28671.900 

-5760.9329 

40547.892 


a For this trajectory, ^ = 70°, i R = 75°, h^ D = 262 n. mi., ty R = 20 n. mi., T = 200 n. mi. , ^ = 0, and k R = 1. 

Plane to centric position vectors are measured in nautical miles. Heliocentric position vectors are measured in astronomical units. 

"Planetocentric vectors are in the planetocentric equatorial (Earth equatorial at Earth, Mars equatorial at Mars) system. Heliocentric vectors are in the heliocentric ecliptic 
system. r 
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TABLE m. - MARS-VENUS-EARTH TRAJECTORY WITH PARKING ORBIT (PERIOD = 80 HR) ABOUT VENUS 
FOR A MARS DEPARTURE DATE OF DECEMBER 15, 1977 (JULIAN DATE 2 443 493) a 


Location 
| (coordinate 

Time from 
departure 


Position vector 
(b, c) 



Velocity vector, fps 
(0 


system) 

periapsis, days 

X 

Y 

z 

Magnitude 

X 

Y 

Z 

Magnitude 

Mars periapsis 
(planetocentric) 

0. 0000000 

-1610.4360 

-487.51092 

-1069.8950 

1993.9524 

1013.3574 

-19876.006 

7531.4149 

21279. 209 

Exit Mars sphere 
of influence 
(planetocentric) 

1.5161525 

107898. 95 

-241735. 97 

164924. 76 

311895.24 

5038.6730 

-10987.367 

7599. 3071 

14277.953 

Exit Mars sphere 
of influence 
(heliocentric) 

1.5161525 

-0.41311593 

1.5511198 

0.04336610 

1.6057765 

-60281.586 

-9563.6723 

3677.5070 

61146. 197 

Enter Venus sphere 
of influence 
(heliocentric) 

259. 57448 

0. 31693437 

-0.64791667 

-0. 02717924 

0.72179090 

121655.43 

58789. 120 

-.69779800 

135295.58 

Enter Venus sphere 
of influence 
(planetocentric) 

259.57448 

-312425.87 

-113050.52 

-24215.915 

333131.74 

19111. 048 

6980.4972 

1082.4660 

20374.769 

Venus periapsis 
(planetocentric) 

260.68412 

1867.3431 

1144.4505 

-2722.0318 

3493.7364 

29696. 826 

7289.7737 

2343.7253 

38527.226 

In parking orbit 

260.68412 

1867.3431 

1144.4505 

-2722.0318 

3493. 7364 

24867.625 

6104.3343 

1962.5964 

32262. 053 

Out of parking orbit 

270.68412 

1867.3088 

1144.4522 

-2722.0546 

3493.7364 

24867.625 

6104.3343 

1962.5964 

32262.053 

Venus periapsis 
(planetocentric) 

270.68412 

1867. 3088 

1144.4522 

-2722. 0546 

3493.7364 

29029.799 

7126.2610 

22910.365 

37661.619 

Exit Venus sphere 
of influence 
(planetocentric) 

271.88744 

99215.993 

-14784. 617 

317670. 21 

333131. 75 

5212.7837 

-977.51682 

17917.953 

18686.404 

Exit Venus sphere 
of influence 
(heliocentric) 

271.88744 

. 52094398 

-.49653745 

-0.03319247 

72044000 

i 

84352.754 

90040. 293 

13425. 145 

124108.32 

Enter Earth sphere 
of influence 
(heliocentric) 

437.82797 

-. 91598292 

.38231432 

. 00423224 

0. 99257587 

-30201. 160 

-84328.989 

-10552. 387 

90193. 354 

Enter Earth sphere 
of influence 
(planetocentric) 

, 437.82797 

-300574.20 

-324349.88 

231981.80 

499362.84 

9383.5571 

9774.8658 

-7263.3681 

15373.863 

Earth periapsis 
(planetocentric) 

440.00000 

2672. 9480 

-141.85234 

-2245.6217 

3493.9361 

2396. 8676 

39336.023 

368. 18051 

39410.698 

a For this trajectory, ^ = 40°, 

iR* 40 ’- \,D = 

150 n. mi. , li 

7T, rl 

= 50 n. mi. , h^ T 

= 200 n. mi. , kj. 

y = 0, and k R = 0. 





b Planetocentric position vectors are measured in nautical miles. Heliocentric position vectors are measured in astronomical units. 


C Planetocentric vectors are in the planetocentric equatorial (Earth equatorial at Earth and Venus, Mars equatorial at Mars) system. Heliocentric vectors are in the heliocentric 
ecliptic system. 
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